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^ ■ ABSTRACT 

^ ; 

, Context. SWIFT J1626. 6-5156 is an X-ray pulsar that was discovered in December 2005 during an X-ray outburst. Although the 

■ X-ray data suggest that the system is a high-mass X-ray binary, very little information exists on the nature of the optical counterpart. 
Aims. We investigate the emission properties of the optical counterpart in the optical and near-IR bands and the long-term X-ray 

I I ' variability of the system in order to determine unambiguously the nature of this X-ray pulsar. 

, Methods. We have performed an X/optical/JR analysis of SWIFT J 1626. 6-5 156. We have analysed all RXTE observations since its 

' discovery, archived optical spectroscopic and photometric data and obtained for the first time near-IR spectra. X-ray energy spectra 

HH ' were fitted with models composed by a combination of photoelectric absorption, a power law with high-energy exponential cutoff and 

(— I ' a Gaussian line profile at 6.5 keV and an absorption edge at around 9 keV. X-ray power spectra were fitted with Lorentzian profiles. 

We identified and measured the equivalent width and relative intensity of the spectral features in the optical and infrared spectra to 

I . determine the spectral type of the optical counterpart. 

O ' Results. The K-band spectrum shows Hel /120581A and HI /121660 A (Brackett-gamma) in emission, which confine the spectral 

\ type of the companion to be earlier than B2.5. The H-band spectrum exhibits the HI Br-18-1 1 recombination series in emission. The 

^ . most prominent feature of the optical band spectrum is the strong emission of the Balmer line Ho-. The 4000-5000 A spectrum contains 

, ^ , ' He II and numerous He I lines in absorption, indicating an early B-type star. The source shows three consecutive stages characterised 
by different types of variability in the X-ray band: a smooth decay after the peak of a large outburst, large- amplitude flaring variability 

, (reminiscent of type I oytbursts) and quiescence. We observed that the spectrum becomes softer as the flux decreases and that this is 

^ . a common characteristic of the X-ray emission for all observing epochs. An emission line feature at ~6.5 keV is also always present. 

I ' Conclusions. The X-ray/optical/IR continuum and spectral features are typical of an accreting X-ray pulsar with an early-type donor. 

' The long-term X-ray variability exhibiting large outbursts, minor outbursts and quiescent emission is also characteristic of hard X-ray 

, transients. We conclude that SWIFT J1626.6-5156 is a Be/X-ray binary with a BOVe companion located at a distance of ~10 kpc. 

\l , Key words. X-rays: binaries - stars: neutron - stars: binaries close 

\6 ■ 
o : 

^ 1. Introduction The fir st detailed study of SWIFT J1626.6-5156 was car- 

ried out bv lRei^ et a l.' (2008), who found that the duration of the 
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SWIFT J1626.6-5 156 is a hard X-ray transient whose actual na- flares was a few hundred seconds and that the X-ray intensity in- 
ture IS uncertain. It was discovered on 18 December 2005 by creased by a factor of 3.5 during the flares. They also presented 



^ Swift/BAT when Its X-ray emission was showing short-term flax- evidence (although not conclusive) in favour of a high-mass X- 
H , ing episodes (|PalmeretaU |2005|). Immediate after the first de- ray binary classification of the system. The determination of the 
. ^ , tection by Swift, the source was observed by RXTE, which con- nature of the system is very important because the flares seen 
firmed the presence of an X-ray pulsar with P^pin = 15.377 jn sWIFT J1626.6-5156 would constitute the shortest events of 
s and strong variations of the pulse fraction during the flares this kind ever reported in a high-mass X-ray binary. 
dMarkwardt & Swankl 120051; iBelloni et al.L 12006). In a subse- 
quent observation with SwiftpiKT, the X-ray position was re- In a recent paper, iBavkal et alJ (l2010l) (see also llcdem et"an 

fined to RA (J2000): 16:26:36.24, DEC (J2000): -51:56:33.5, '2011') performed a pulse frequency analysis and reported the 
with a 90% error radius of 3.5" (Campana et al., 2006). The discovery of the orbital parameters of the system. They also 
optical and infrared observations of the only counterpart so far studied the evolution of the neutron star spin period through- 
proposed a re contradictory. While optica l spectroscopy suggests out the outburst. The periodic trend of pulse frequencies yielded 
a Be star dNegueruela & Marcol l2006h. infrared o bservations an orbital period of PqA - 132.89 + 0.03 days and an ec- 
seem to indicate a late-type object (" Reaet all 12006). The y-ray centricity of e = 0.08 + 0.01. Pulse-phase spectroscopy has 
mission INTEGRAL also detected SWIFT J1626.6-5156 about also been addressed bv iReig etal . (2008), while a preliminary 
a year after i ts dis covery, but no X-ray flares were observed analysis of the evo lution of the pulse profiles can be found in 
dTarana et alil2006l) . iBavkal et al.l (l201 ih . Here we present a detailed X-ray spectral 
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Table 1. Optical/infrared apparent magnitudes of the optical 

counterpart to SWIFT J1626.6-5156 from the NOMAD cata- Fig.2. H- and K-band spectra of the infrared counterpart to 
logue. Note that B, V, R and / are photographic magnitudes. SWIFT J1626.6-5156. Helium and hydrogen lines are cleai'ly 

seen in emission. 



B 


V 


R 


I J 


H 


K 


16.81 


15.54 


15.81 


14.27 13.44 


12.95 


12.54 



and timing analysis of the entire 2006 outburst and subsequent 
flaring emission. The X-ray varibility that followed the major 
outburst is reminiscent of the so-call ed type-I outburst, seen in 
many Be/X-ray binaries (see e.g. Wil son et al.Ll2002h . However, 
unlike typical Be/X-ray binaries, the frequency of the outbursts 
observed in S WIFT J 1626.6-5 15 6 does not always equal the 
orbital period dBavkal et al.Ll2010b . For this reason, we shall re- 
fer to these quasiperiodic enhancements of the X-ray intensity 
as flares ( not to be confused with the three short-term flares re- 
ported bv'Re ig et al.l ([2008')). We also analysed data during the 
X-ray quiescent state that followed the flaring state and present 
new optical and infrared observations. Our aim is to solve the 
nature of this unusual accreting X-ray pulsar. 



Table 2. H and A'-band line identifications and measured equiv- 
alent widths for 2MASS J16263652-5156305. 



Spectral 
feature 


Wavelength 
dim) 


Equivalent 
width (A) 


Br- 10 


1.737 


-7.8+0.3 


Br- 11 


1.678 


-7.9+0.5 


Br- 12 


1.641 


-8.3±0.5 


Br- 13 


1.611 


-7.0+0.5 


Br- 14 


1.588 


-9.7+0.6 


Br- 15 


1.570 


-8.3+0.5 


Br- 16 


1.556 


-7.8+0.2 


Br- 17 


1.544 


-8.0+0.3 


Br- 18 


1.535 


-6.2+0.3 


He I 


2.058 


-9.9+0.4 


Bry 


2.166 


-8.8+0.3 



2. Optical and infrared observations 

2.1. NIR Observations 

Figure [1] shows a digital image of the field around the X-ray 
best position with the identification of the optical counterpart. 
Near-IR data were obtained in visiting mode in July 2010, at 
the European Southern Observator y (ESO). The employed in- 
strument was the SofI spectrograph dMoorwood et alii 19981) . on 
the 3.5m New Technology Telescope (NTT) at La Silla, Chile. 
We used the long slit spectroscopy mode, at low resolution 
(R = 588) with the 1.53-2.52/imgrismand 1" width slit. The in- 
strument large field objective provided a FOV of 4.92' X 4.92'. 
The sky had thin cirri and seeing averaged between 0.9" and 
1.3". 

With the aim of ensuring accurate removal of atmospheric 
features from t he spectrum, w e follo wed a strategy similar to 
that oudined bv lClark & Steele! (12000). At the telescope, we ob- 
served an AOIV standard star immediately before the target and 
a G3 V immediately after it, in order to obtain very small diff'er- 
ences in airmass (differences between 0.02 and 0.04 airmasses 
were accomplished). To compute the telluric features in the re- 
gion of the HI 21 661 A (Brackett-y line, or Bry), which is the 
only non-telluric feature in the A-star spectra, we employed the 
observed G-star spectra divided by the solar spectrurrQ properly 
degraded in resolution. The dispersion solution obtained for the 
SofI spectra was also applied and the spectra of the A star, G 
star and the solar one were aligned in wavelength space. For the 
K spectral region, a telluric spectrum was obtained by patching 
into the A-star spectrum the ratio between the G star and the 
solar spectrum in the Bry region (we selected the range 21 590 - 
21 739 A). For the H spectral region, we employed as the telluric 
spectrum, the ratio between the G-star and the solar spectrum. 



' We used the NSO/Kitt Peak FTS solar spectrum, produced by 
NSF/NOAO. 
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Table 3. Log of the optical spectroscopic observations 



Date 


MJD 


instrument 
set-up 


Grism 
number 


Wavelength 
range(A) 


EW(H 
(A) 


a) 


EW(H/3) 
(A) 


15-02-2006 


53782.359 


EMMI/NTT 


grism#6 


5700-8700 


-40 ± 


2 




17-09-2007 


54361.008 


EFOSC/3P6 


grism#18 


4700-6770 


-45 ± 


2 


-5.0 ±0.9 


17-09-2007 


54360.984 


EFOSC/3P6 


grism#14 


3095-5085 






-4.3 ± 0.5 


18-09-2007 


54362.025 


EFOSC/3P6 


grism#14 


3095-5085 






-4.2 + 0.5" 


18-09-2007 


54361.993 


EFOSC/3P6 


grism#3 


3050-6100 






-4.5 ± 0.5* 


01-06-2009 


54983.078 


EFOSC/NTT 


grism#13 


3685-9315 


-45 ± 


2 


-5.1+0.9 



a: Average of three measurements 
b: Average of two measurements 



Data reduction was performed using the IRAF0 package, fol- 
lowing the standard procedure. We first corrected for the inter- 
quadrant row cross-talk, a feature that affects the SofI detector; 
we then applied sky subtraction; we employed dome flat-fields 
and extracted the one dimensional spectrum. Wavelength cali- 
bration was accomplished using Xenon and Neon lamp spectra. 
Spurious features, such as cosmic rays or bad pixels, were re- 
moved by interpolation, when necessary. The reduced spectra, 
one covering the H band, one the K band, were normalized by 
dividing them by a fitted polynomial continuum. We finally cor- 
rected for telluric absorption, dividing each scientific spectrum 
by its corresponding telluric spectrum, obtained as described 
above. A scale and a shift factor were applied to the telluric spec- 
trum, to best correct for the airmass difference and the possible 
wavelength shift; the optimum values for these parameters were 
obtained using an iterative procedure that minimizes the residual 
noise. The final H and K spectra are shown in Fig.|2l 

2.2. Optical spectra 

Optical spectroscopic data were retrieved from the ESO archive 
facility. Table [3] shows the log of the optical spectroscopic ob- 
servations, summarising the instrumental set-up and giving the 
equivalent width of the Ha and HyS Unes. The EMMI obs e rvatio n 
corresponds to that reported by iNegueruela & Marco! (l2006h . 
while the EFOSC observations are unpublished: programs IDs 
079.D-0371(A) and 083.D-0110(A). Data reduction was per- 
formed using the IRAF packages, following the standard pro- 
cedure. 



3. X-ray observations and data analysis 

We a nalysed data obtained by all three instruments aboard RXTE 
(Uradt et al., 1993): the All sky Monitor (ASM) data consist 
of daily flux averages in the energy range 1.3-12.1 keV. The 
Proportional Counter Array (PCA) covers the lower part of the 
energy range 2-60 keV, and consists of five identical coaligned 
gas-filled proportional units giving a total collecting area of 6500 
cm"^ and provides an energy resolution of 18% at 6 keV. The 
High Energy Timing Experiment (HEXTE) is constituted by 2 
clusters of 4 Nal/Csl scintillation counters, with a total collect- 
ing area of 2 x 800 cm^, sensitive in the 15-250 keV band with 
a nominal energy resolution of 15% at 60 keV. Data taken dur- 

^ IRAF is distributed by the National Optical Astronomy 
Observatories which is operated by the Association of Universities for 
Research in Astronomy, Inc. under contract with the National Science 
Foundation. 



ing satellite slews, passage through the South Atlantic Anomaly 
and Earth occultation were removed. The overall on-source time 
analysed amounts to 496.3 ks. 

The mission-specific packages of Heasarc FTOOLS (version 
6.6.3) were employed to perform data reduction, while XSPEC 
vl2.43 was used for spectral analysis. For each observation we 
obtained an average energy spectrum and a power spectrum. 

The energy spectra were generated from standard mode data, 
namely, the PCA Standard! of VCm and HEXTE "FS58" clus- 
ter B data, which have a time resolution of 16 s and cover the 
2-60 keV range with 129 channels and the 15-250 keV with 64 
channels, respectively. All spectra were background-subtracted. 
A systematic error of 0.6% was added in quadrature to the PCA 
spectra to account for systematic errors. The power spectra were 
generated by computing Fourier transforms of 128-s segments 
of the 2"^-s binned light curves in the energy range 2-15 keV 
(PCA channels 0-35) and averaging them together The final 
power spectra were logarithmically rebinned in frequency and 
corrected for dead time effects according to the prescription 
given in Nowaketal. (1999). Power spectra were normalized 
such that the integral over a given frequency range equals to 
the squar ed fractional r ms amplitude jBelloni & Hasingenl 19901 ; 
Mivamot o"et al.Lll99lh . 

To fit the energy spectra, we used a model composed by a 
combination of photoelectric absorption, a power law with high- 
energy exponential cutoff, a Gaussian line profile at 6.5 keV and 
an absorption edge at around 9 keV to account for Fe K fluo- 
rescence. We did not find evidence for a cyclotron absorption 
feature, commonly seen in other accreting X-ray pulsars. The 
power spectra were fitted using Lorentzian profiles only. The 
main peak of the X-ray pulsations and two to three of its har- 
monics are clearly seen in the power spectra. These peaks were 
fitted with narrow Lorentzians (the width was normally fixed at 
~ 0.001 Hz). The broad-band noise was fitted with two zero- 
centred Lorentzians. The brighter observations required one ex- 
tra Lorentzian to account for the noise above ~ 7 Hz. 

4. The optical counterpart 

4.1. Spectral type 

The information on the optical counterpart to SWIFT J1626.6- 
5156 is very scarce. To the authors' knowledge no optical pho- 
tometric observations are published. The only dedicated ob- 
servations are tho se reported in the Astronomer's Telegrams 
ATel#713 (iRea et a l.. 2006) who obtained infrared phot omet- 
ric magnitudes and ATel#739 dNegueruela & Marcoll2006li) . who 

^ http://heasarc.gsfc.nasa.gov/docs/xanadu/xspec/ 
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Fig. 3. Optical spectrum in the classification region. The identified lines are He I /14009, /14026, /14144, /i4387, /i4471, /14713 A; 
He II ^4541 and ^4686 A; Si IV i4089 A; the C III+OII blend ^4650 A. 



performed spectroscopic observations and reported the presence 
of Ha strongly in emission. 

Information on photometric magnitudes and colours can 
be found in astronomical catalogues. We searched for pos- 
sible optical/infrared counterparts in the NOMAE Q catalogue 
around the best-fit X-ray position (ICam pana et al.L |2006). The 
optical counterpart to SWIFT J1626. 6-5156 is the star named 
as NOMAD 1-0380-0705569 (=2MASS 16263652-5 156305 = 
USNO-Bl.O 0380-0649488) and it is shown in Fig. □ The 
selected s t ar is the only one exhibiting significant i nfrared 
jRea et al.L l2006h and Ha (iNegueruela & Marcd l2006l) emis- 
sion within the SWIFTPCRT error circle or its vicinity, which are 
the signatures of the vast majority of optical counterparts in ac- 
creting X-ray pulsars. Table [T]gives the catalogued photometric 
magnitudes. 

The low-resolution /T-band spectrum is almost featureless, 
except for the presence of He I /120581 A, slightly affected by 
a residual telluric absorption and HI /(21 660 A (Brackett-y, or 
Bry), which appear strongly in emission. These two lines are 
typical of Be-stars and confine t he spectral type of the compan- 
ion star to be earlier than B2.5 ( ClcU-k & StedS. |2000). The H- 
band spectrum is more complex, with the H I Br- 18-1 1 recombi- 
nation series in emission. As for the K band, all the observed fea- 
tures are typical of Be stars, while they are absent, or observed 



The NOMAD catalogue is a facility provided by the U. S. Naval 
Observatory and contains astrometric and photometric data for about 
1.1 billion stars derived from the Hipparcos, Tycho-2, UCAC2, Yellow- 
Blue 6, and USNO-B catalogs for astrometry and optical photometry, 
supplemented by 2MASS near-infrared photometry. The primary aim 
of NOMAD is to help users retrieve the best currently available astro- 
metric data for any star in the sky by providing these data in one place. 



in absorption in normal (non-emission) OB stars. The observed 
features, with measured equivalent widths are shown in Table|2l 

By comparison with spectral atlases (" Clark & Steelel 120001 : 
Steele & Clark, 2001), the optical counterpart to SWIFT 
J1626. 6-5156 can be classified as a B0-2Ve star, confirming the 
nature of the system as a high-mass X-ray binary. Note that all 
spectroscopically identified optical companions of Be/X-ray bi- 
naries in the Milky Way have spectral types in the narrow range 
09-B2, with a peak around BO (see e.g. Fig. 6 in Antoniou et aU, 
I2l)09l) . 

Based on JHK photometrv.' Rea et alJ 12006') suggested that 
the optical counterpart to SWIFT J1626.6-5156 could not be an 
early-type star because after de-reddening using the X-ray hy- 
drogen column density (A^h = 0.9 x 10^^ cm"^), the infrared 
magnitudes and colours were too faint to be consistent with 
an OB star. The value of A^h correspo nded to a Swift/XRT ob- 
servations reported bv ICampana et al.l (.2006) . How ever, the ac- 
tual value of A^H reported in ICampana et al. (l2006l) is an order 
of magnitude larger than that used bv lRea et all (l2006h . We re- 
peated the computation of the intrinsic c olour of the source, em - 
ploying the SWIFT/XRT Nh value from lCampana et alJ (l2006l) . 
A^H = (9.4 + 1.0) X 10^' cm"^, which is more reliable than our 
RXTE measurement due to the softer energy range sensitivity. 
This implies a visual interstellar extinction Ay = 4.25 + 0.45 
mag (Giiver & Ozel, 2009), which converts into an interstellar 
IR color excess E'' ( H - K ) = 0.063 xAy ^ 0.27 ± 0.03 mag 
(iRieke & Lebofskvl Il985h . Note that the dense circumstellar 
material surrounding the Be star introduces extra r eddening that 
cannot be ignored, especially at long wavelengths (iDachs et al.L 
1988; Dougherty et al., 1994). Hence E""{H - K) = E''(H - 
K) + E'^^{H - K). The reddening caused by the disc can be es- 
timated from the relationship between infrared excess and the 
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Fig. 4. X-ray light curve of SWIFT J1626.6-5156 showing five 
years worth of data. Each data point corresponds to the average 
count rate of the individual pointings. Also shown are the soft 
colour, SC=4-7 keV/2-4keV and hard colour, HC=10-15 keV/7- 
10 keV. 



Ha equivalent width from iHowells et al.l (1200 ll) . who derived 
E^'iH - K) ^ Qm6EW{Ha) - 0.030 mag. We do not know 
the state of the disc when the infrared observations were mea- 
sured, but if we assume that the Ha equivalent width was at the 
level reported here ~ -40 A, then E"{H - /T) = 0.2 1 + 0.02 mag 
and the resulting intrinsic colour is {H - K)q - {H - A')2mass - 
E'°'{ H -K) = -(0.07 ± 0.04) mag, consistent with an early-type 
star (lKoornneeflll983l) . The errors were obtained by propagating 
the errors of the measurements. 

The September 2007 blue-end spectra allow us to refine the 
spectral type estimated from the near-IR data. Figure[3]shows the 
optical spectrum of SWIFT J1626. 6-5156 in the traditional clas- 
sification region (3800-4800 A). Note that emission affects the 
lines of the Balmer series: HyS and Hy are in emission, while WS 
is partially filled-in with emission. The spectrum contains many 
He I fines (.14009, /14026, /14144, .14387, .14471, A4713 A), in- 
dicating an early B-type star. The presence of Hell (.14541 and 
.14686 A) indicates that the spectral type of the optical star is ear- 
lier than B 1 . Specifically He II A4686 A is last seen at B0.5-B0.7 
dmibom & Fitzpatrickl|l993). On the other hand, the weakness 
of He II .14541 A relative to He I i4471 A indicates an spectral 
type later than 09. 

With regard to the luminosity classification, the weakness of 
Si III .14552-68 A and the fact that the line intensity of He II 
/14686 A is larger than that of He 1 .14713 A seem to indicate a 
main-sequence star Other luminosity indicators are the ratios of 
He lines to nearby metallic lines. SWIFT J1626. 6-5156 shows 
He II A4686/CIII .14650 > 1 and He 1 .14144/Si IV 4089 ~ 1 . An 
evolved star would have those ratios < 1 and < 1 , respectively. 

A visual comparison of the SWIFT J1626.6-5156 spec- 
trum (Fig. [3]l with those of MK standards in the atlas by 
(IWalborn & Fi tzpatrickl 1 19901) reveals that the spectrum resem- 
bles that of the standard star v Ori, a BOV star. Therefore, we 
conclude that the optical counterpart to the X-ray accreting pul- 
sar SWIFT J1626.6-5156 is a BOVe star. 




10 100 
PCU2 count rate (c/s) (2-15 keV) 



Fig. 5. Hardness-intensity diagram. 



4.2. Distance 



Although the uncertainty of the spectral type may introduce 
some extra uncertainty in the distance estimation, the main 
source of error in estimating the distance stems from the scat- 
tering of the intrinsic colours and absolute magnitudes associ- 
ated to each spectral and luminosity class of the calibrations. 
Even for the same spect ral type, varia tions of up to 1 .5 mag in 
the absolute magnitude (I Weened l2006i) and of 0.04 mag in the 
intrin sic (B - V)q colours are found i n studies from various au- 
thors (|JohnsonL ll966VWegneii [T994l) .lJaschek & Go med (Il998h 
analysed the absolute magnitude of about 100 MK standards and 
concluded that the intrinsic dispersion of the mean absolute mag- 
nitude amounts to 0.7 mag. We have adopted this value as the 
error on the absolute magnitude. 

The BVRI magnitudes listed in Table [1] are photographic 
magnitudes, extracted from the USNO-Bl and the YB6 sur- 
veys, and hence they are not directly comparable with the intrin- 
sic colours and absolute magnitude calibrations of the Johnson- 
Cousins photometric system. This prevent us from using these 
magnitudes to compute the interstellar reddening and the dis- 
tance to our target. We will approach this task by using the 
2MASS JHK magnitudes also listed in Table [1] 

Assuming that the optical counterpart to SWIFT J1626.6- 
5156 is a BOV star, then (J - K)o - -0.17 mag (Koornnee| 
Il983h . The total IR excess is E""(J - K) ^ 1.07. Part of this ex- 
cess is due to the circumstellar continuum emission of the Be 
star. We have calculated the circumstellar contribution to the 
IR excess to Be E'^^{ J - K) - 0.25 mag from the equivalent 
width of the Ha line and the H- and K- band spectral features , 
by means of the formulae presented by iHowells et al.l (1200 ih . 
Thus we obtained the interstellar IR excess as £"(7 - K) - 
E""iJ - K) - E'HJ - K) = 0.82 mag. TTiis irn pfies a visual 
extinction of Av = 4.7 mag (iKoornneefl Il983h . little higher, 
but consistent with the value obtained from the X-ray observa- 
tions. TTi£_absolutemagnitude of a BOV star is My - -3.95 
mag (IWegneit '2006), which can be converted into Mk since 
{V - K)o - -0.85 mag (Koornneef, 1983); the intrinsic mag- 
nitude Ko was computed as the difference of the 2MASS mag- 
nitude and the K-band extinction A/^ resulti ng from the rel ation 
AJE(J -K) ^ 2.4(.l)-i ", for A = 2.2^m (iDraind [T989I) . We 
finally calculated the distance from Mk - Kq + 5 - 5logd, ob- 
taining d = 10.7 + 3.5 kpc. This value must be considered as a 
lower limit, since we neglected the circumstellar contribution to 
Mk. The error was estimated assuming an error in the observed 
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infrared magnitudes (hence E{J - K)) of 0.03 mag and 0.7 mag 
in the absolute magnitude. 

5. Long-term X-ray variability 

The long-term X-ray hght curve of SWIFT J1626.6-5156, 
covering five years worth of data, since its discovery on 19 
December 2005 (MJD 53723) up to 14 December 2010 (MJD 
55544), is shown in Fig.|4] This figure also shows the evolution 
of the X-ray colours defined by the ratio of the count rates in the 
energy range 4-7 keV over 2-4 keV (soft colour, SC) and 10-15 
keV over 7-10 keV (hard colour, HC). The X-ray flux of SWIFT 
J1626. 6-5156 decreased exponentially during the first 90 days 
of the observations. On 20 March 2006 (MJD 53814) the source 
began to exhibit low-amplitude X-ray flares whose duration var- 
ied between 3-8 days. Nevertheless the overall decrease in flux 
continued until 26 June 2006 (MJD 53912), where a large flare 
took place. After the large flare, the mean intensity increased to 
~10 c s ' PCU ' and the X-ray flaring behaviour became more 
apparent. A total of 15 flares can be seen in Fig|4] The time dif- 
ference between the peak of the flares varied from ~45 days of 
the first half to -9 5 days of the second half (Reig et al., 2008; 
iBavkal et al.Ll2010l) . 

On 29 November 2008 (MJD 54799) the X-ray count rate 
decreased suddenly by almost an order of magnitude to ~ 1 .5 c 
s ' PCU ' and the large-amplitude flaring pattern ceased. The 
source entered a quiescent phase. Note that the source was ob- 
served daily until MJD 53827, every 2-4 days until MJD 53936, 
every 9-10 days until MJD 54900 and in a more irregular way 
every 15-20 days (but also including a few closer observations) 
until the end of the observations on MJD 55544. 

The soft colour follows closely the evolution of the X-ray 
flux (Fig. Is). As the count rate increases so does the SC. The 
HC also shows a positive correlation with flux. However, this is 
weaker than that at lower energies and it is evident only at higher 
count rates >20 c s ' PCU ' . The correlation seems to extend to 
the lower count rate data but the scattering of the data points due 
to poor statistics is too large below ~ 3 c s"' PCU"' . 

Given the richness in X-ray variability exhibited by SWIFT 
J 1626.6-5 156, we have defined three different time intervals and 
performed a separated analysis on each one. These periods cor- 
respond to the smooth decay part of the outburst, including the 
small-scale flaring period (MJD 53723-53912), the large scale 
flaring (MJD 53912-54799) and the quiescent interval (MJD 
54799-55544). 

5.1. Outburst decay (MJD 53723-53912) 

SWIFT J1626.6-5156 was discovered on 18 December 2005 
when the source was near the peak of a major outburst. The X- 
ray flux of the first observations was 4.5 x lO"** erg cm"^ s"' and 
decreased for the next ~190 days (interval A in Fig.|4|i. Figure|6] 
shows two representative energy and power spectra of the out- 
burst decay at fluxes diflFering by about one order of magnitude. 
The spectral parameters changed smoothly as the outburst de- 
cayed (Fig.|7]i. However, at flux below Fbreak ~ 1.5 x 10"^ erg 
cm"^ s"' the source experienced a significant change in its spec- 
tral shape. The X-ray continuum of the brighter observations, 
i.e., when the 3-30 keV flux was above Fbreak, showed a faster 
decay at high energies requiring an exponential cutoff to fit the 
data. The energy of the cutoff increased as the flux decreased. 
Below Fbreak this Component was not needed and the X-ray con- 
tinuum was well represented by a single absorbed power law 




Energy (keV) Frequency (Hz) 



Fig. 6. Energy and power spectra at high (upper panels) and low 
(lower panels) X-ray flux. The X-ray flux is ~ 4 x 10"^ and 
~ 5 X lO^'^erg cm"^ s"', respectively. 

extending to energies up to 100 keV (Fig.|6]l. Likewise, the hy- 
drogen column density, A^h fell rapidly as the flux decreased be- 
low Fbreak, whilc it remained at a constant value of ~ 6.5 x 10^^ 
cm"^ above f break- The photon index decreased as the flux in- 
creased, that is, when the source was bright the X-ray emission 
was harder. At low flux, roughly coincident with Fbieak the pho- 
ton index flattened at a value of 1 .7 (Fig. |7]l. Near the peak of the 
outburst the photon index was -1.1. 

The presence of the iron emission line or the absorption edge 
was not always statistically significant. When the line parameters 
could not be well constrained, the line and/or edge energy were 
fixed at their average values, prior to the fit. The mean value of 
the central and absorption edge energies were calculated using 
observations in which the intensity of the line (normalization) 
was more than 2cr significant and resulted in =6.5 +0. 1 keV 
and /iedee =9.3 +0.4 keV, where the errors represent the standard 
deviation of all the measurements. The three bottom panels of 
Fig. Q show the dependence of the iron line energy and intensity 
and the edge energy with the 3-30 keV X-ray continuum flux. 
There is a tight correlation between the X-ray continuum and 
the intensity of the iron line, indicating that as the illumination 
of the cool matter responsible for the line emission increases, so 
does the strength of the line. 

The characteristic frequency of the broad-band noise in- 
creased as the flux increased, but saturated at v w 1 .5 Hz above 
^^break- The large error bars at low flux weakens the statistical 
significance of this trend. Almost identical trend was followed 
by the fractional amplitude of variability. Above ~ 2.6 x 10"^ 
erg cm"- s"' there is excess power at high frequencies and an 
extra noise component with characteristic frequency at ~ 6 - 7 
Hz and nns 10-15% was needed to obtain acceptable fits (upper 
right panel in Fig.|6]l. Some of the brightest observations contain 
a weak QPO with centroid frequency at ~ 1 Hz and nns ~ 3%. 
Although the inclussion of this component improves the fit, an 
F-test shows that the probabiUty that this improvement occurs 
by chance is ~l-3%. 

5.2. Large-scale flaring emission (MJD 53912-54799) 

At the end of the outburst decay and after a series of 
low-amplitude flares, SWIFT J1626. 6-5156 began to display 
quasiperiodic increases in the X-ray intensity (interval B in 
Fig©- The maxima and minima of these flares remained roughly 
at the same flux level, with Fmax = (28 + 5) x 10"" erg cm"^ 
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Fig. 7. Evolution of the spectral parameters over the decay of the outburst. Flux corresponds to the 3-30 keV energy range. Power- 
law index, hydrogen column density and cutoff energy (top panels) and central energy, edge energy and intensity of the iron line 
(bottom panels) as a function of X-ray flux. 

Table 4. Average best-fit spectral parameters during the large-amplitude flares and the quiescent state that followed. 



Parameter 



Minimum 



Flaring state Quiescent state 

rise maximum decay I 



Photon index (F) 
Cutoff energy* (£cuc) 
Normalization** 
Line energy* (£fc) 
Line width* (crpe) 
Edge energy* (£edge) 
Edge optical depth (Tcdgc) 



1.9 + 0.1 

30-' 

0.017 ± 0.001 
6.9 + 0.1 
0.4^ 

9-' 

07+°"' 

"•"'-0.07 

1.9/48 



1.5 + 0.1 
20+5 
0.031 +0.003 
6.7 + 0.1 

0.3^ 
8.5 + 0.2 
0.17 + 0.05 
1.6/46 



1.16 + 0.09 
11 + 1 
0.034 + 0.003 
6.6 + 0.1 

0.2^ 
8.3 + 0.1 
0.26 + 0.03 
1.1/46 



1.6 + 0.1 
21+^ 

0.032 + 0.004 

6.7 + 0.1 
0.3^ 

8.5+0.3 
0.17 + 0.06 
1.2/45 



2.10 + 0.08 

0.007 + 0.001 
6.6 + 0.1 
0.1^ 



1.0/36 



Notes. *: keV; **: ph keV 'cm ^ s ' at 1 keV; ^ : fixed 



s"' and Fmin = (5 + 2) x 10"" erg cm"^ s"'. The only excep- 
tion was the first flare, which started and reached ~2.5 times 
lower and higher flux than the average minimum and maximum 
flux, respectively. The timescale of the X-ray modulation var- 
ied between 45 and 95 days, although most of the time sepa- 
ration between flares gathered around Po rb/2 and f orb/3, being 
f orb — 133 days (see lBavkal et al.L 1201 Ol for an analysis of this 
behaviour). 

Due to the low count rate, the energy spectra above ~10 
keV and the power spectra above ~ 1 Hz are too noisy for a 
meaningful analysis. The cutoff energy and especially the iron 
emission line parameters cannot be constrained in the individual 
spectra, which represent data obtained by integrating over 0.5-1 
ks, typically. However, by joining observations with roughly the 
same flux, those components become significant. Therefore in 
order to increase the signal-to-noise ratio we obtained an aver- 



aged peaked, trough, intermediate-rise and intermediate-decay 
flux energy and power spectrum. The total on-source time for 
these average spectra was 17.6, 19.1, 15.4 and 11.5 ks, respec- 
tively. Table |4] gives the results of the spectral fits. The spec- 
trum at the peaks is harder than that of the troughs. The cut- 
off energy follows the same trend with X-ray flux as that seen 
during the outburst decay, namely, it increases as the flux de- 
creases. The characteristic frequency of the broad-band noise 
component also agrees with the values that the source showed 
during the decay of the outburst for the same flux level. At the 
peak of the flares v^ax = 0.40 ± 0.09 Hz, while at intermediate 
flux Viise ~ Vdecay = 0.12 + 0.05 Hz. The fractional amplitude of 
variability is ~ 25% in all cases. The timing parameters at the 
minimum flux of the flares cannot be constrained due to the low 
signal-to-noise ratio of the power spectrum. 
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Fig. 8. Poib -^'spin and Poii,-EW(Ha) diagrams. The star symbol 
marks the position of SWIFT J1626. 6-5156. Open symbols in 
the top panel correspond to the peculiar SGXBs 2S 0114+65 
and OAO 1657-41 and the BeXB SAX J2 103 .4+4545. 



5.3. Quiescence (MJD 54799-55544) 

After the series of flares the source entered a quiescent state, 
characterised by little variability. The X-ray flux decreased to 
an average level of (1.7 + 0.4) x 10"" erg cm"- s"'. Again, 
the individual spectra are too noisy to be able to constrain the 
spectral parameters. Neither an iron emission line nor the expo- 
nential cutoff are statistically required in the individual spectra. 
However, when an average energy spectrum is obtained by com- 
bining all observations in the interval MJD 54799-55544 then a 
good fit cannot be obtained unless a Gaussian component at 6.5 
keV is included, although no edge is required. The inclusion of 
that component reduced the statistic from 94 for 38 degrees 
of freedom to 35 for 36 degrees of freedom. An F-test gives a 
probability that the improvement of the fit occurs by chance of 
only 1.9 X 10"^. The total on-source time of this spectrum was 
46.4 ks, corresponding to 61 observations. The best-fit spectral 
parameters are given in Table |4] The energy range was limited to 
3-20 keV. The spectrum during the quiescent state is softer than 
during the outburst decay and during the flaring state, but fol- 
lows the general trend also observed in the other two states that 
the lower the flux, the larger the photon index, i.e., the softer the 
spectrum. 



6. Discussion 

The available data in various wavelength bands provide strong 
evidence that the system is a Be/X-ray binary. In the X-ray band, 
the detection of 15.4-s pulsations and the shape of the spectral 
continuum (a power law with exponential cutoff) are characteris- 



tic of accreting X-ray pulsars. In the optical band, the strong Ha 
emission line and the presence of He I and He II lines indicate 
an OB-type star as the optical counterpart. In the NIR, both the 
H- and the A'-band observed features (Hel end HI recombination 
series in emission) are also typical of early Be stars. This allows 
us to classify the counterpart to SWIFT J 1626.6-5 156 as a BOVe 
star, and the system as a Be/X-ray binary. 

Further evidence in favour of a Be/X-ray bin ary cl assifica- 
tion comes from the f oi-b - Psmn dia gram (Corb e3.[T986l) and the 
Pa,b-EW{Ha) diagram (iReig et al.LTr997:,ReigLl201 ih . FigurelSl 



shows these two diagrams. As can be seen, SWIFT J1626.6- 
5156 (big star symbol) agrees with the f m-b - EW{Ha) correla- 
tion and falls in the Be/X-ray binaries region of the Porb - ^'spin 
diagram. 

The optical spectra that cover the region of the Ha line show 
comparable values of the Ha equivalent width despite the fact 
that they correspond to very different X-ray states and were taken 
years apart. The February 2006 observation corresponds to the 
middle of the X-ray outburst decay, the September 2007 obser- 
vations correspond to the middle of the flaring episode and the 
June 2009 observation to the beginning of the quiescent X-ray 
state. In Be stars, the strength and shape of the Ho- line provides 
information about the physical conditions in the circumstellar 
disc. Generally, after a major X-ray outburst, the Ha equivalent 
width decreases substantially indicating a weakening of the disc 
(Coe et al., 1994) or even the line shape turns from an emission 
profile to a n absorption profile, suggesting that the disc is lost 
(iReig et all 12007,) . In contrast, the disc in SWIFT J1626.6-5156 
does not seem to be affected by the X-ray outburst, as the Ha 
equivalent width did not vary much immediately after the X-ray 
outburst nor a few years later 

The behaviour of SWIFT J1626.6-5156 resembles that of 
the Be/X-ray binary KS 1947+300. As SWIFT J1626.6-5156, 
KS 1947+300 went through a series of periodic increases of 
the X-ray intensity (t ype I outburst) after a major (type II) out- 
burst (see Fig. 11 in lReigll201 1.) without the Ha line noticing 
it. Another characteristic in common with KS 1947+300 is the 
near-circular orbit. The eccentricity in these two systems (0.03 
and 0.08, respectively) are among the smallest in Be/X-ray bina- 
ries. An alternative scenario would be that oflA 1118-61, which 
exhibi ted a huge equiva lent width (-90 A) before the X-ray out- 
burst dCoe et all Il994l) . After the outburst, the Ha equivalent 
width decreased substantially but remained at a relatively high 
level around -60 A for the following years (Villada etal., 1999*). 

A discrete component at ~6.5 keV is always present, albeit 
weak, in the X-ray energy spectra of SWIFT J 1626.6-5 156 in all 
states. This component is interpreted as reprocessing of the hard 
X-ray continuum in relatively cool matter Near neutral iron gen- 
erates a line centered at 6.4 keV, while this energy increases as 
the ionisation stage increases. However, the energy separation is 
so small that even ionised iron up to Fe XVIII can be thought as 
part of the 6.4 keV blend (Liedahl, 2005). The line energy during 
the outburst decay is consistent with this neutral iron blend. The 
line energy during the flaring state seems to be larger. However, 
given the weakness of the line we do not claim a higher ioni- 
sation degree in this state. What it is important to notice here 
is that the fluorescent iron line provides strong evidence for the 
presence of material in the vicinity of the X-ray source. It has 
been detected in virtually all high-mass X-ray binaries with su- 
pergiant companions. In these systems the Fe line fluorescence is 
produced in the stellar wind of the massive star (iTorrejon et all 
2010). The possible site for reprocessed emission in a Be/X-ray 
binary are an accretion disc or the circumstellar (decretion) disc. 
The former surrounds the neutron star, while the latter surrounds 
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the Be star's equator. The detection of this component in the 
energy spectrum of the quiescent state would favour the equa- 
torial disc as the site for the reprocessing of hard-energy pho- 
tons, as we would expect that the accretion disc, if ever present, 
would have disappeared. The alternative scenario that the Fe iron 
line emission is produced in the thin thermal h ot plasma that is 
presumably located along the Galactic plane ( Yamauch i et al.L 
[2009) is ruled out by the observed correlation between the line 
intensity and the X-ray continuum flux. 

We estimate the distance to be ~ 10 kpc with an error of the 
order of 30%. This distance and the Galactic coordinates of the 
source (I - 332.8, b - -2.0), place the system in the Norma Arm 
of the Galaxy. The sources of error in the determination of the 
distance can be attributed to: i) uncertainties in the calibrations 
that give intrinsic colours and absolute magnitudes as a function 
of spectral types, ii) the presence of a large circumstellar disc 
which contaminates the optical and IR colours making the ob- 
served colours appear redder that their intrinsic values, and to a 
lesser extent. Hi) the uncertainty in the spectral type. 

Note that in Be stars the total reddening not only includes 
the contribution of the interstellar extinction but also of the 
reddening due to the circumstellar material. The contribution 
of the circumstellar disc to the photometric colours increases 
with wavelength from ~0.1 ma g in (B - V) to ~ 0.8 mag in 
(J - M) for well-developed discs (Da chs et aUll988h . Moreover, 
it has been shown that the circumstellar reddenin g correlates 
with the equivalent width of the Ha lin e (Dachs et al.L Il988t 
iFabregat & Torrei6nLll998tlHowells et al.Li2001.) . The larger the 
equivalent width, the larger the contribution of the circumstel- 
lar emission to the total reddening. Since SWIFT J1626.6-5156 
shows strong Ha emission with an equivalent width of ^5 A, 
the disc emission can make an important contribution to the 
colour excess. Although we t ry to a ccount for this extra ex- 
tinction using the equations by lHowe lls et al. (2001), it is worth 
noticing that the photometric infrared magnitudes and the opti- 
cal and infrared spectra were not contemporaneous, which intro- 
duces extra uncertainty in the distance estimation given the high- 
amplitude variability that characterises the optical and infrared 
emission of Be stars. Our dist ance estimat e is con sistent with the 
value of ~15 kpc reported bv llcdem" et all (1201 lb . These authors 
estimated the distance fro m the correlation b etween spin-up rate 
and X-ray flux. Although llcdem et al.l (1201 ih do not give the er- 
ror of their measurement, the application of this method in other 
sources has shown that the uncertainty is of the order of ~10- 
30%. 

7. Conclusion 

We have solved the nature of the X-ray pulsar SWIFT J1626.6- 
5 156. The new optical and infrared data provide strong evidence 
that the system is a Be/X-ray binary. The location of SWIFT 
J1626. 6-5156 in the PqA - ^'spin and forb - EW{Ha) diagrams, 
the X-ray pulsations, the shape of the X-ray continuum, the long- 
term X-ray variability, the emission of the Balmer and Brackett 
lines in the optical/IR spectrum are typical characteristics of 
high-mass X-ray binaries with a Be companion. We conclude 
that SWIFT J1626.6-5156 is a Be/X-ray binary with a BOVe op- 
tical counterpart located at a distance of ~ 10 kpc. Optical photo- 
metric observations are needed to reduce the uncertainty in the 
distance estimation. 
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